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ABSTRACT. We describe here the Instrumentation, data, and prelimin- 
ary results from the summer, 1984, solar oscillation observing 
program which was carried out using the 60-foot tower telescope of 
the Mt. Wilson Observatory. This program was carried out with a 
dedicated solar oscillation observing system and obtained full-disk 
Dopplergrams every 40 seconds for up to 11 hours per day. Between 
June and September, 1984, observations were obtained with a Na 
magneto-optical filter on 90 different days. The data analysis has 
progressed to the point that spherical harmonic filter functions have 
been employed to generate a few one-dimensional power spectra from a 
single day*s observations. 

1. INTRODUCTION 

Shortly after the initial observations of high-degree solar p- 
mode oscillations (Deubner, 1975; Rhodes, Ulrich, and Simon, 1977; 
Deubner, Ulrich, and Rhodes, 1979) two of us (EJR and RKU) recognized 
the need for a dedicated solar oscillation observing system. 
Accordingly, we began to develop such a system Jointly with Dr. Robert 
P. Howard at the Mt. Wilson Observatory. The Initial design of this 
system was developed during 1978 and major modifications were made to 
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thim design faring 1980 cad 1981. The two ant significant las>rove- 
•cuts were the replacement of ooe'dlaeasloul let icon detectors with 
* CID csaera (jstea, sad the woe of a magneto-optical filter Instead 
of a spectrograph. The developmental states of this system was 
first described by Khodes et al. (1981). More recent reports about 
this system were given by Ihodes et al. (1983 and 1984a). The first 
solar oscillation poser spectrum obtained with this systen during 
July, 1983, was presented in Ihodes al. (1984a and 1984b). 

Bere we wish to report on the flr&. extended operation of this 
oscillation system in the manner for which It was orglnally designed: 
the acquisition of daily full-disk Dopplergraas with high time 
resolution, faring the fall of 1983 and the spring of 1984 the full- 
disk Imaging capability of the system was Installed and tested with 
the goal of obtaining such observations throughout the 1984 observing 
season. He will first report on the observations which were obtained 
and then finish with a progress report on the status of the reduction 
of that data. The analysis of a larger quantity of this data is 
currently In progress with our primary goals being an laproveaent In 
the published accuracy of low- and intermediate-degree p-mode frequen- 
cies and the probing of solar Internal rotation as functions of depth 
and latitude as Brown has reported earlier at this meeting (Brown, 
1985). 

2. OBSERVATIONS 

The observations obtained during the summer of 1984 consisted of 
pairs of narrow- band filtergrams which were obtained in the two Na 
D lines. The filtergrams were obtained with a two-cell version of 
the magneto- opt leal filter (MOF) which was originally Invented by 
Cacdanl. This filter was first described In Clmino et al. (1968a 
and 1968b), and In Caccianl et al. (1970). Its uses for solar ob- 
servations were later described in Agnelli, Caccianl, and Fofl 
(1975); In Caccianl and Fofl (1978); in Caccianl, Fortini, and 
Torelll (1980); and in Caccianl et al. (1981). More recently, its 
principles of operation were reviewed by Caccianl and Rhodes (1984), 
who described both one-and two-cell versions of It. 

Figure 1 shows the two- cell configuration of the MOF which was 
employed during May through September, 1984, to obtain the data 
described here. During that time the filter employed two 1000-Gauss 
permanent magnets. (Subsequent to September, 1984, the first of 
these two magnets was replaced with a 6000-Gauss permanent magnet and 
this combination of 6000- and 1000-Gauss magnets has been adopted for 
all subsequent Na observations.) 

At the left of the magnets In Figure 1 we show the mechanism 
which controlled the moveable quarter-wave plate. This mechanism 
consisted of a reversible stepper motor and a Geneva movement. When 
operated under the control of an LSI-11 microcomputer, this mechanism 
allowed us to rotate the wave plate through angles of +45* and -45° 
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Figure 1. Magneto-Optical Filter. ThiB figure shows the two-cell Na MOF an it was Installed at Mt 
Wilson during 1984. At that time the filter employed two 1000-Gauss magnets. 



with respect to the crossed polarizers which were slso part of the 
WF. 

When the save plate was positioned in the +45* position, for 
example, we were able to record filtergrams in the bandpass located 
on the red side of the HaD lines, while a rotation of the wave plate 
to the -45* angle tuned the filter to the blue bandpass* 

The exposures of the Photometries, ltd. CID cawera which we 
employed were regulated by a rotating disk shutter* The rotation of 
this shutter was governed by a series of uniformly-spaced stepper 
aotor pulses which were supplied to the shutter by a Superior Electron- 
ic Buffered Translator* Each flltergraa took 5 seconds for both 
exposure and cawera read-out and so a pair of red- and blue-wing 
flltergraas could have been obtained once every 10 seconds. In fact, 
this sequence was only repeated once every 40 seconds for up to 11*4 
hours per day. Therefore, a single day's observing run consisted of 
a waviiwM of 1024 red flltergraas and 1024 blue flltergraas. 

The transmission profiles of the two - cell HOF were similar to 
those shown in Figure 2. Here we show two more-recent prfiles which 
were obtained in the NaDl line with the 6000-6au8s/1000-Gau8s combina- 
tion of magnets mentioned earlier. These wavelength scans were 
obtained with the spectrograph of the 150-foot tower telescope 
at Mt. Wilson. The two profiles corresponding to the +45* and -45* 
positions of the quarter-wave plate are superimposed upon each other 
in this Figure. The degree of separation of the two bandpasses is 
Illustrated by the low level of the transmission of each curve at the 
location of the peak of the other curve. 

The cumulative amount of data obtained during our 1984 campaign 
is shown in Figure 3. Here the number of hours of cloud-free 
observations obtained on each day is plotted as a function of time. 
The initial observations were obtained on May 8, while the last 
observations were obtained on September 9. The stlpled portion of 
the diagram shows that most of the May data was subsequently destroyed 
when all but the first one of each day's tapes were later re-written 
with newer data in early August. A single tape per day from that 
period was saved in order to have some Dopplergrams for eventual 
comparison with those produced at the 150-foot tower on the same 
days. On some of the days between May 20 and June 8 several other 
tests were run with the M0F and with the observing hardware; hence, 
the absence of data during that interval was not due to a period of 
bad weather. 

Figure 3 shows that on several occasions during the campaign we 
did Indeed obtain 7 to 11 hours of cloud-free data per day for 
intervals exceeding one week in duration. On only 12 days after with 
June 24 was the weather bad enough to prevent the recording of any 
data at all. The fact that this much data was obtained during 1984 
was due solely to the diligence of S. Tomczyk who observed and archived 
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the data almost single-handedly for the entire three months. 

A histogram of the number of days having a given length run is 
shown in Figure 4. Inspection of this figure shows that on 45 of the 
90 days at least 8 hours of data was obtained. Since this histogram 
only Includes the number of cloud-free hours on each day, the total 
amount of data that was obtained actually exceeded the amount shown 
In the histogram. The large number of days having between one and 
two hours of data was the result of our effort to record at least a 
single tape's worth of data even when It became obvious early In a 
given day that the weather would not permit the full -day run to be 
completed. 

3. CONVERSION INTO DOPPLERGRAMS 

At the conclusion of each day's observations the time series of 
alternating red and blue filtergrams was archived from disk to tape. 
The first steps in the analysis of this data were the computation of 
the location of the disk center of each filtergram and the transla- 
tional registration of each filtergram to a common center location. 
The registration process was necessary to correct for the Image 
motion which occurred between the acquisition of each red filtergram 
and its corresponding blue companion. It also allowed us to correct 
for 8 low drifts in the solar image introduced by the gulder on the 
time scale of several hours. 

After the filtergrams were registered in the computer, each pair 
could then be converted into a Dopplergram with the following simple 
formula: 


*blue(*»y) ” *red(*»y) 

velocity (x,y) - constant* (I) 

*blue( x »y) ^ *red(*»y) 

where Ibiue an< * *red referred to the intensities of the pixel at lo- 
cation (x,y) in the blue and red filtergrams and the constant was a 
calibration constant which had to be determined in some way. 

A Dopplergram similar to those employed with the 244x248 pixel 
CID camera and the Na MOF is shown here in Figure 5. The figure 
shows a Dopplergram which was obtained with a JPL 800x800 pixel CCD 
camera (although only an area 500 pixels on a side was employed in 
this instance). This Dopplergram was obtained on June 1, 1984. It 
was obtained by differencing two filtergrams which were taken 15 
seconds apart; hence, the appearance of the solar limb is due to 
slight image distortions caused by the Earth's atmosphere during and 
between the exposure of the two filtergrams. Solar rotation is 
visible as the shading from white to black across the visible disk. 

By comparison we show the Dopplergram obtained in X5250 with 
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Figure 4. Mean Acquisition Histogram. The number of days having a given run length are plotted as 
function of observing run length. 
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Figure 5. MOF Dopplergram. Na Dopplergram obtained with MOF and 
800x800 pixel CCD camera on 1 June 1984. Solar rotation is visible 
as the shading across the disk. 



the spectrograph of the 150-foot tower on the same day in Figure 6. 

Again, solar rotation Is visible scross the disk as are several 
active regions. The velocity flows around the active regions are less 
pronounced In the A5250 DoppLergram than they were In Figure 5. 
Some of this difference may be due to the fact that the NaD lines are 
chromospheric lines while X5250 ij a photospherlc line. Additional 
tests at Mt. Wilson will be needed to clarify this possibility. 

The constant given in equation (1) can in principle be deter- 
mined in two different ways and hence confirmed Independently. 
First, it is possible to use the rate of solar rotation as measured 
in a particular latitude zone ci a given day at the 150-foot tower to 
provide the radial velocities against which the observed intensity 
differences can be compared. The idea is illustrated schematically 
in Figure 7. Here the scaled Doppler signal (defined for now as 
50000 • [Ij - Ig]/[Ig + I R J) from a full-disk Dopplergram obtained 
with the CID camera is shown as a function of position along a swath 
that was located parallel to the solar equator. If the MOF were re- 
sponding linearly to the solar radial velocity over a range of ± 200 
m/sec, then we would expect to see a perfectly straight line running 
from one end of the equator to the other, while the slope of this 
straight line would be related to the solar rotational velocity. 
That is, for the 150-foot tower Dopplergram of Figure b, a slope of 
roughly 2 m/sec/arcsecond would be obtained near the equator. 

A regression analysis of the so-called Doppler signal against 
the measured solar velocities along the equator would then yield a 
straight line whose slope could later be used to convert a particular 
intensity difference anywhere on the disk into an equivalent radial 
velocity. The degree of constancy of this calibration constant over 
the solar image would then be a measure of the presence or absence of 
systematic velocity errors introduced by the MOF. The principal 
objection to employing this calibration method is that it makes use 
of the 150-foot tower Dopplergrams to provide measurements of the 
very quantity which the MOF is trying to measure, namely, the solar 
radial velocity at each point on the visible disk. 

An alternative calibration procedure is to compute an average 
Doppler signal for every single Dopplergram and to then compute what 
the predicted Earth-sun radial velocity was for the moment in time. 
This later calculation can be made with an ephemeris program and 
subsequently a regression analysis can be carried out of these Doppler 
signals as a function of the predicted velocity. The slope and 
Intercept of the regression line that is so obtained will then be the 
calibration constants necessary for conversion of the intensity 
differences into radial velocities. 

Two different applications of this latter approach are illustra- 
ted here in Figures 8 and 9. In the upper left panel of Figure 8 we 
have plotted the mean Doppler signal (10000 [Ig-Ig * IIb+IrJ) which 



Figure 6. Comparison Dopplergram. A Dopplergram obtained In A5250 
with the spectrograph at the Mt . Wilson 150' tower telescope on 1 
June 1984. 




we obtained by averaging the raw Doppler signal fro* individual 
pixels within a 512^ x 512^ area at the center of the solar disk. 
Here this wean Doppler signal is shown as a function of the Universal 
tine at which each frane of data was recorded. 

In the upper right panel of Figure 8 we have plotted the radial 
velocity of the Earth relative to the sun as a function of Universal 
Tine. The radial velocities shown in the panel were coaputed with an 
enpheneris program which was originally developed at JPL for the 
Interplanetary program, and which includes the effects of planetary 
perturbations on the Earth. The output from this JPL program was 
combined with computations of the Earth’s diurnal motion and of the 
Earth's lunar barycentric motion in order to provide the total radial 
velocity between Mt. Wilson and the sun at every moment of time. The 
rotation of the Earth is clearly vlcible In this panel. 

In the bottom panel we have plotted the mean Doppler signal versus 
the corresponding computed sun-Earth velocity for that moment of 
time. A linear regression analysis was performed upon the data shown 
in this panel and we obtained the regression coefficients listed at 
the top of the Figure. The slope of 1.222 m/sec per Doppler unit Is 
what would later be used to convert individual pixel difference 
signals into radial velocities. The Intercept of the regression line 
Includes the effects of: 1) the gravitational redshift, 2) asym- 
metries in the solar line profile, and 3) asymmetries in the trans- 
mission profile of the MOF, as were seen In Figure 2. The correla- 
tion coefficient of 0.993 from the regression analysis shows that the 
Doppler/velocity conversion was Indeed quite linear for this particu- 
lar data. 

The MOF data used to generate Figure 8 were obtained In July 
1983. When we repeated this analysis with the 1984 full-disk data 
described in Section 2, we got the results shown in Figure 9. Here 
the mean Doppler signal shown in the upper left panel was obtained by 
averaging the Doppler signal over the entire visible disk of the sun. 
Again the computed sun-Earth radial velocities for the observing 
times of each frame are plotted at the upper right. And again, at 
the bottom, we show the results of plotting the mean Doppler signal 
against the computed radial velocity. Here the relationship is not 
nearly as linear as was that shown in Figure 8. While we are not yet 
certain of the cause of this departure from linearity, we have noticed 
that the discrepancy is worst at large hour angles early in the 
morning and late in the afternoon. Thus, we currently believe that 
telluric H 2 O contamination in the wings of the NaD2 line, which would 
have the greatest effect at the largest air masses, may be introducing 
the time-dependent departure from linearity which we are seeing in 
Figure 9. We will repeat this type of analysis for every day’s worth 
of data In order to see whether or not this departure from linearity 
is always present in the 1984 data. (In the meantime, we have observed 
only in the D2 line during 1985 since that line has much less telluric 
contamination. ) 
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Figure 9, Re-Callbratlon Using Full-Disk Averages. The calibration procedure shown In Figure 8 Is 
repeated here with Doppler signals obtained from full-disk averages of Dopplergrams taken on 27 July 
1984. Here the departure from linearity Is believed to be due to time-varying telluric contamination 
of the Na D2 line. 




4. INITIAL DATA ANALYSIS 


4.1 Spherical Harmonic Filtering 

The next step In the data reduction sequence was the spherical 
harmonic filtering of the various time series of Dopplergrams. In 
order to carry out this filtering operation, we developed a computer 
program with which we could project spherical harmonic functions onto 
the planar 240 x 240 pixel grid of the CID camera Images. Since the 
image registration program which translated all of the Dopplergrams 
from a given run to a common center also computed the radius of each 
solar image In pixels, we needed a program which was general enough 
to correct for the varying size of the solar image. In addition, we 
wanted a program which would allow for the position angle of the 
solar rotation axis and for the non-zero heliographic latitude of the 
sub-earth point on the solar disk. We also felt that it would be 
desirable for the harmonic generation scheme to be accurate up to 
1-values of 1000 or more. Thus, the program which we developed can 
now generate fore-shortened spherical harmonics of the appropriate 
size and of any orientation. 

Once the spherical harmonic generator had been developed, we 
next used its output to perform the spatial filtering. This was 
accomplished by dot-multiplying each Dopplergram with a given spheri- 
cal harmonic on a plxel-by-pixel basis. This dot-multiplication 
corresponded to a double integration of the product of the Doppler- 
gram and the harmonic. The output of this double integration was a 
separate filtered time series for each individual harmonic. 

An example of such a filtered time series is shown here in Figure 
10. In this Figure we show the detrended amplitude of the £ ■ 10, m 
■ 10 sectoral spherical harmonic as a function of time on July 27, 
1984. Because of the non-linearity of the time-dependent calibration 
analysis tor tnis July 27 observing run, as was shown in Figure 9, we 
did not in fact subtract off the computed sun-Earth velocity from 
each Individual filtered Dopplergram. Instead, we decided to employ 
a simple third-order polynomial detrending analysis to remove the 
diurnal motion of the Earth from the filtered time series until we 
could correctly determine how to calibrate the July 27 dataset. 
Therefore, the time series displayed in Figure 10 has been both 
filtered and detrended. 

Once the detrending analysis was complete, the information which 
remained in the filtered time series shown in Figure 10 was the 
amplitude which resulted from the beating together of all of the 
various solar oscillation modes having the spatial pattern given by 
the £ - 10 sectoral harmonic, as well as some cross-talk from other, 
nearby sectoral and near-sectoral harmonics. (An analysis is now 
underway to determine more precisely the amount of cross-talk expected 
from such adjacent spherical harmonics.) 





4.2 Power Spectral Computations 

Once a filtered time series had been generated it was then a 
simple matter to employ a one-dimensional FFT program and co^ute the 
corresponding one-dimensional power spectrum. Three such power spec- 
tra from the July 27, 1984 run are shown in Figures 11, 12, and 13. 
Figure 11 shows the power spectrum which resulted from transforming 
the filtered time series shown in Figure 10. Therefore, the relative 
power levels of several 1 » 10 sectoral harmonics (along with some 
cross-talk from adjacent harmonics) are Illustrated in this Figure. 
The power levels shown on the vertical axis are in arbitrary units 
since we had not employed the correct nc malization constants at the 
time this spectrum was originally compute-.. 

In Figure 11 the frequency unit employed is the circular fre- 
quency * “ Instead of the linear frequency, v. In these units the 
peak of the "5-minute" p-mode band is near w • 0.020 sec' 1 . Due to 
the limited frequency resolution available from a single 11-hour 
observing run, the individual p-modea are not fully resolved in this 
spectrum. Nevertheless, the relatively low noise level of the MOF 
can be seen at frequencies between 0.004 and 0.016 sec 1 and above 
0.060 sec' 1 . 

At this time we should point out that the filters were operated 
without the benefit of any temperature stabilization scheme whatso- 
ever. Furthermore, they were operated within the spectrograph cham- 
ber without the presence of any insulating material to shield them 
from ambient air temperature fluctuations. Lastly, the power sup- 
plies which controlled the heaters on the MOF cells, and hence the 
optical depth of the Na clouds within the cells, were not operated 
with any power regulation beyond that built into the supplies them- 
selves. The rise in power level at frequencies below w ■ 0.004 
sec' 1 could be due in part to uncontrolled drifts in the power inputs 
to the cell heaters. (A closed-loop servo system which will control 
the temperature of the Na pits within the MOF cells to ± 0.1 °C i a now 
nearing completion and will soon be installed to test these ideas.) 

Two more power spectral slices from the same day’s data are 
illustrated here in Figures 12 and 13. In Figure 12 the power spectrum 
wa3 completed from the dataset that had been filtered with the A ■ m 
- 20 sectoral harmonic, while in Figure 13 the A • m « 40 was used. 
The power levels are again not expressed in (m/sec) 2 , although both 
power spectra should be normalized consistently with the spectrum of 
Figure 11. (That is to say that, while the absolute value of the 
normalization coefficients may not be correct, their relative values 
are. All of the factors which are A* and m-dependent have been 
correctly Included in the normalization coefficients already.) Note 
that while the noise level looks higher in Figures 12 and 13 than in 
Figure 11, the vertical scales are different. In fact, the absolute 
noise levels are quite similar in all three cases. 
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Figure 13. Sane as Figures 11 and 12 except that this spectrun corresponds 
sectoral harmonic. 
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TASKS REMAINING 


The tasks which remain In our processing of the summer 1984 date 
can be devided Into three major areas. First, there are several 
difficulties with the prelieinary processing of the data which we 
have already outlined above and which we must resolve before we can 
proceed. Second, there Is the nonaous task of actually carrying out 
the analysis on the large number of data tapes which we have accumu- 
lated. Third, there Is the scientific interpretation of the frequen- 
cies and frequency splittings which we will obtain from the power 
spectra. 

There are four principal difficulties which comprise the first 
group of tasks. These are: 1) the problem with the non-linearity of 
the time-dependent velocity calibration analysis, 2) the search for 
possible spatial non-unlformltles in the calibration coefficients, 3) 
the determination of the correct spherical harmonic normalization 
scheme, and 4) the re-coding of the spherical harmonic generation 
program and the spatial filtering program for our CSPI MAP 300 and 
Mini-MAP array processors In order to cut their execution times down 
to a manageable level. 

The actual production runs to be carried out in the second phase 
will Include: 1) the registration and calibration of the Individual 
Dopplergrams for many more days of observations, 2) the spatial 
filtering analysis, 3) the time-synchronization of the filtered time 
series for multi-day runs, 4) the application of gap-filling tech- 
niques for the nighttime data gaps, 5) the confutation of the power 
spectra from these multi-day filtered time series, and 6) the deter- 
mination of the modal frequencies from the resulting power spectra. 

The key scientific analyses to be carried out include: 1) the 
confutation of the absolute frequencies for all of the modes observed 
as a function of harmonic degree, 2) the comparison of these frequen- 
cies with the predictions of various solar models, 3) the computation 
of the modal frequency splittings Introduced by solar rotation, 4) 
the inversion of the frequency splittings in order to measure inter- 
nal solar rotation as a function of both depth and latitude, and 5) a 
comparison of the absolute frequencies and the frequency splittings 
with those from earlier epochs In order to search for possible time- 
dependent changes In the sun. 

6. SUMMARY 

He have presented a review of the acquisition and preliminary 
analysis of the most extensive set of imaging data yet obtained 
anywhere for the study of solar oscillations. He have outlined the 
unique aspects of the instrumentation which have made the acquisi- 
tion of this much data possible. He have gone on to describe some 
currently-unresolved difficulties with the data (as well as methods 




which will hopefully Allow us to resolve these difficulties). And 
we have presented e few preliminary results fro* our analytes of only 
a small amount of the available data. Finally, we have outlined the 
tasks which remain to be completed in our analysis* 
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